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Abstract 

The projected properties of triaxial generalization of the modified Hubble mass 
models are studied. These models are constructed by adding the additional ra- 
dial functions, each multiplied by a low-order spherical harmonic, to the models 
of Chakraborty & Thakur (2000). The projected surface density of mass models 
can be calculated analytically which allows us to derive the analytic expressions of 
axial ratio and position angles of major axis of constant density elliptical contours 
at asymptotic radii. The models are more general than those studied earlier in the 
sense that the inclusions of additional terms in density distribution, allows one to 
produce varieties of the radial profile of axial ratio and position angle, in particu- 
lar, their small scale variations at inner radii. Strong correlations are found to exist 
between the observed axial ratio evaluated at 0.25i? e and at AR e which occupy well- 
separated regions in the parameter space for different choices of the intrinsic axial 
ratios. These correlations can be exploited to predict the intrinsic shape of the mass 
model, independent of the viewing angles. Using Bayesian statistics, the result of a 
test case launched for an estimation of the shape of a model galaxy is found to be 
satisfactory. 
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1 Introduction 

Observed photometric properties of majority of elliptical galaxies show isopho- 
tal twists and variations in eliipticity and position angle of major axis with 
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radius. These observed propertie s can b e produced in a natural way by triaxial 
models. Chakrabortv fc Thakur ( 2000L hereafter CTOO) studied the projected 



properties of a family of mass models which are triaxial g eneralisation of modi- 
fied Hubble mass model. The model was first proposed bv lSchwarzschild ( 19791 ) 



as a numerical model for a triaxial st ellar system in dynamical e quilibrium and 
later casted into analytical form by de Zeeuw fc Merritt ( 1983| ). It was found 



by CTOO that these models show ellipticity variations and isophotal twists in 
their projections along the line of sight. Moreover, the radial profiles of the 
parameters of the elliptical isophotes were found to be smooth functions of 
the radius (see CTOO). However, many elliptical galaxies, devoid of any fea- 
tures, indicating the absence of shells or dust, are found to exhibit small scale 
variations in the radial profiles of the parameters of elliptical isophotes. This 
indicates that density distributions of such ellipticals may be more complex 
than that of CTOO. 

Above fact inspired us to modify the models of CTOO by adding the addi- 
tional radial functions, each multiplied by a low-order spherical harmonic. 
The projected surface density £ of resultant mass model can be calculated 
analytically. This makes possible to investigate some of the projected proper- 
ties analytically. We calculated the profiles of surface density S, axial ratio b/a 
and position angle G* of the major axis as functions of radius. The analytical 
expression of the £ is very useful for the calculations of b/a and 9*. In the 
asymptotic regions analytical expressions can be derived, while in the interme- 
diate regions simple numerical methods can be adopted for these calculations. 
The inclusions of additional terms in the density distribution, allows us to 
produce varieties of the radial profile of b/a and B* that can be compared 
with photometric data of real galaxies. 



It was shown by de Zeeuw fc Carollol ( 199d ) that the observed ellipticities and 



the position angle, at small and at large radii can be used to derive the intrinsic 
axial ratios of the mass model, as a function of the viewing angles. Since the 
viewing angles of galaxies are largely unknown, it will be worth-while to get 
an estimate of the intrinsi c axia l ratios independent of the viewing angles. In 



this re gard. iStatler fc Frvl (|1994 hereafter SF94), and lThakur fc Chakrabortv 
( 200lL hereafter TC01) found that in case the observed parameters exhibit cor- 
relations when a model with a given set of intrinsic parameters is viewed in 
all possible orientations, one can obtain a probable estimate of the intrinsic 
shape, independent of the viewing angles. Furthermore, it was also shown by 
TC01 that the intrinsic shapes of triaxial mass models can be estimated us- 
ing photometric properties and further, the results of shape estimation would 
be insensitive to the choice of models, if one considers ensembles of models 
which represent above mentioned correlations between the observed parame- 
ters. Thus, another aim of our present study is to build ensembles of models 
showing the correlations between observed parameters so that one could get 
an estimate of model independent shape. 
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In §2, we describe the mass model and in §3, we present the projected prop- 
erties. Our results are described in §4, and §5 is devoted to summary and 
discussion. 



2 The mass model 



We modified the triaxial potential, given in de Zeeuw fc Merritt ( 1983| ). corre- 



sponding to the density distribution considered in CTOO by adding the addi- 
tional radial functions v e (r) and w e (r), each multiplied by a low-order spherical 
harmonic. Thus, our assumed potential V(r, 9, 0) has a following form 

V(r, 9, 0) = u(r) - (v(r) + v e (r))Y 2 °(9) + (w(r) + w e (r))Y 2 2 (6, 0), (1) 

where u(r), v(r), w(r), v e (r) and w e (r) are five radial functions, (r, 0) are 
spherical coordinates defined such that x = r sin 9 cos 0, y = r sin 9 sin and 
z = r cos 9, and Y 2 {9) = ^cos 2 9 — ~ and Y 2 2 (9, 0) = 3sin 2 9cos2(f) are the usual 
spherical harmonics. 

We take u{r) to be the potential of the spherical modified Hubble mass model, 
defined by 



ln\r + Jb 2 n + r 2 ] 

u{r) = -GM— (2) 

r 

where M is the total mass of the model and b n is the scale len g th. W e choose 
v(r) and w(r) to be same as adopted in de Zeeuw fc Merritt (jl983h and re- 



produced below with slightly different notations. 
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where b\, ...,64 are constants. The additional radial functions v e (r) and w e {r) 
are chosen such that they are effective in intermediate range of r only and not 
at small and at large radii. In order to satisfy this condition, a suitable choice 
of v P (r) and wJr) are 



GM afr 6 



bl (a 2 + r 2 ) 4 
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, . GM afr" 

where ai, 04 are constants. We note that v(r) and iu(r) go as — r 2 at small 
radii and as — 1/r, at large radii. On the other hand, t> e (r) and u> e (r) go as 
— r 6 at small radii and as —1/r 2 at large radii. 

For the above considered potential, the associated density distribution p(r, 9, 0) 
follows from Poisson's equation 



p(r, 9, 0) = f(r) - (g(r) + g e {r)) Y 2 °(9) + (h(r) + h e {r)) Y 2 2 (9, 0). (5) 
The five radial functions appeared in equation (5) are given by 
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The radial dependence of additional functions g e (r) and h e (r) in the density 
distribution are such that at asymptotic radii, the density remains same as 
that of CTOO. At large radii, g e (r) and h e (r) decrease as r -4 , whereas g(r) 
and h(r) decrease as r~ 3 . Likewise, at small radii, g e (r) and h e (r) decrease as 
r 4 , whereas g(r) and h(r) goes as r 2 . Thus, the radial functions f(r), g(r) and 
h(r) are dominating at small and at large radii in the density distribution (5). 
The four ratios (bi/b Q ), ... , (6 4 /6 D ), appeared in g(r) and h(r), are expressed 
in terms of axial ratios of the approximately ellipsoidal constant p surfaces 
at small and at large radii by the same way as adopted in CTOO. Prescribing 
the values of axial ratios of the approximately ellipsoidal constant p surfaces 
at very large and at very small radii as (poo^oo) an d (p ,Qo), respectively, we 
find that 
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Here it is clear that the expressions of (bi/b a ) and {b^/b ) are exactly similar 
to those of CTOO, whereas they are found to be completely different for (62/^0) 
and (64/60)- Moreover, throughout the paper, the constant parameters a\, 0,4 
are chosen such that |g e (r)| <C \g(r)\ and |/i e (r) | <C \h(r)\ at all r. 



3 Projected properties 



3.1 Projected surface density 



The density form (5) of the models allows a straight forward calculation of the 
pr ojected surface den s ity £ for any viewing angle. We adopted the convention 
of de Zeeuw fc Franxl (jl,989) to project the model along line-of-sight. We chose 
coordinates (x , y , z ) such that z -axis runs along the line-of-sight, and x - 
axis lies in the (x, y) plane. We considered (6 , <p ) as the standard spherical 
coordinates of the line-of-sight and (R, O) as the polar coordinates in (x , y ) 
plane. The projected surface density £(i?, 6) is obtained by integrating the 
model density along the line-of-sight, 



+00 



£(i?,e)= J pdz, (; 

—00 

and after doing some intermediate calculations, it takes the following form 



£ (R, Q) = E (R) + £ 2 (#) cos 2 (6 - 0,; 
where 



(9) 



£ D (R) = 2Fi + 1 - 3 cos 2 9 



(Gi + G\ e ) — ^(G 2 + ^2e 



+ [6(^i + H le ) - 9(# 2 + H 2e )\ sin 2 0' cos : 



S 2 (i?)= 6 (H 2 + H 2e ) cos9' sin: 
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3(# 2 + H 2e ) (1 + cos 2 6) cos 



(10) 
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We have defined the integrals 



G 1 (R) = [ ^ 9{r) dT , Gi e (R) = [ \ 9e{r) dr , 
{ yV 2 - W) I ^(r 2 - R 2 ) 



oo 

G 2 {R) = R 2 j 

R 



g(r) dr 



r ^{r 2 - R 2 ) ' 



G 2e (R) = R 2 J 

R 



9e{r) dr 
r ^(r 2 - R?y 



and similarly integrals for F±, Hi, H 2 , H\ e and H 2e in terms of functions /(r), 
h(r) and h e (r), respectively. 

Above integrals are calculated analytically and can be expressed as 



F 1 (R) 
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h = p(63R 4 + UR 2 a 2 + 3a 4 ), 

/ 2 = (3(Q3R 6 + 21RW + 9/?V + 3a 6 ), 

h = (3(63R 8 + 28R 6 a 2 + 18i? 4 a 4 + 12/? 2 a 6 + 7a 8 ), 

/ 4 = /3(63 J R 2 + 7a 2 ). (13) 

Similarly, Hi(R) and H 2 (R) can be written by substituting 63 and 64 in place 
of bi and 6 2 in G\(R) and G 2 (R), respectively. On the other hand, Hi e (R) and 
H 2e (R) are defined by substituting a 3 and a 4 in place of a\ and a 2 in G\ e (R) 
and G 2e (R), respectively. 

The projected surface density (9) has its major axis at position angle 0*, 
which is given by 
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ta " 2e - ftl + a-W < 14 > 

where hi, h 2 and h 3 depend only on the viewing angles (#', 0') and are defined 

as 



h\ = sin 2 — cos 2 cos 2 6 , 
h 2 = cos 2 — sin 2 4> cos2 $ > 

/i 3 = sin 20 cos# . (15) 

The quantity T = T(R) in equation (14) is the triaxiality parameter, which 
is given by 



T = T(m 4 {H 2 {R) + H 2e {R)) 

~ 1 ] (G 2 (R) + G 2e (R)) + 2 (H 2 (R) + H 2e (R))' 1 J 

Equation (14) gives the position angle 0* of the major axis when it satisfies 
the condition 



(H 2 + H 2e )h 3 sm2e*<0. (17) 

The observed axial ratio b/a of the projected surface density can be calculated 
by using the relation 



E(&,e*-|) = E(a,0*), (18) 
which can be rewritten as 



S (a) + S 2 (a) = S (6)-S 2 (6). (19) 

At a finite radial distance from the centre, the axial ratio can be calculated 
numerically using equation (19). However, at asymptotic limits, the analytical 
expressions can be obtained for b/a. 



3.2 Projected properties at large and at small radii 



The ratio (G 2 + G 2e )/(H 2 + H 2e ) in (16) can be calculated analytically using 
(12), which in turn enables us to write an analytical expressions for T and 0*. 
In particular, at large R, the ratio \{G 2 + G 2e ) / \H 2 + i?2e)]oo and the triaxiality 
parameter are given by 
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(20) 
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(21) 



which clearly indicate that they are similar to {G 2 /H 2 ) 00 and triaxiality pa- 
rameter at large R of CT00, respectively. Thus, the position angle 6* at large 
R would also be identical to that of CT00 and can be calculated using equation 
(21) in (14). 

Likewise, at small R, the ratio [(G 2 + G 2e )/(H 2 + H 2e )} and the triaxiality 
parameter T Q are 



(G 2 + G 
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It is clear that [(G2 + G 2e )/(H 2 + i^2e)]o and T D differ much from (G 2 /H 2 ) 
and triaxiality parameter at small R of CT00, respectively, which implies that 
the position angle 6* at small R would also be much different than that of 
CT00. 

We have used equation (19) to calculate the axial ratios at asymptotic radii. 
In order to write it in some particular form, we define a few more functions 



9 

h = - (h + h 2 ) 2 , 
h 6 = 9 {h\ - hi) 
of viewing angles (#',0 ) and then further define 



(24) 



Z = 



h 4 + /i 5 



G 2 + G 2e 
H 2 + H 2e 



+ h 6 



G 2 + G 2e 
H 2 + H 2e 



1/2 



(25) 



8 



and 
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At very large R, the axial ratio is given by 



b\ 2 _ bl-2bl Zoo 



where is the value of Z at very large R which can be evaluated by substi- 
tuting the value of [(G 2 + G 2e )/(H 2 + # 2e )]oc from (20) in (25), i.e., 



Zoo = 



*•+*• (I) +h ° (I 



1/2 



(28) 



On the other hand, at very small R, it can be expressed as 



C - 12B% 
K_ 

C+12B% 



(29) 



Thus, we find that the position angle 9* of major axis and the axial ratio b/a 
are the same as those of CT00 at large radii, whereas they are much different 
at small radii. The reason behind this fact is that light from all radii will 
contribute to the projected properties evaluated at small radii and hence, the 
effects of additional radial functions are seen there. 
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4 Results 



4-1 Radial profiles of the projected properties 

The projected surface density £ of our considered triaxial mass model (5) can 
be calculated analytically. Its approximate elliptical isodensity contours show 
the variations in the axial ratio b/a and the position angle 6* of the major 
axis as functions of R. These features can be clearly seen from Fig. 1 where we 
present the radial profiles of £, b/a and 6* for particular choices of intrinsic 
axial ratios, p Q = Poo = V and q Q = q^ = g, the viewing angles {6 , ) and the 
constant parameters (04, ). On the basis of these results, our model 

can be considered as an example of simple analytical model showing ellipticity 
variations and isophotal twists in its projection along the line-of-sight. Since 
the resultant density distributing given in (5) can have different contributions 
of the additional radial functions g e and h e for the various choices of constant 
parameters (a±, 02, 03, 04), in Fig. 2 we plotted the profiles of b/a as a function 
of R by allowing the different values of these constant parameters in such 
a manner that the resultant p must remain positive. In this way, the effects 
of additional radial functions on the profiles of b/a can be understood. We 
found that varieties of b/a profile can be produced by changing the constant 
parameters (01,02,03,04) whose trends are found to be different from each 
other. Besides, it has also been noticed that the inclusions of additional radial 
functions produce small scale variations in b/a profiles, specially at inner R, 
as compare to the smooth b/a profiles found in case of CTOO. Moreover, the 
profiles of position angle G* also show small scale variation at inner R for 
different choices of the constant parameters (04,02,03,04) and in Fig. 3, we 
have shown it only for one case which has all the intrinsic parameters similar 
to those considered in bottom right panel of Fig. 2. The isodensity contours of 
many real elliptical galaxies show such small scale variations in the profiles of 
b/a and G* as functions of R rather than the smooth monotonically increasing 
or decreasing profiles of b/a and Q* as those found in CTOO. This suggests that 
the radial profiles of b/a and O* of our present traixial models, as a function 
of viewing angles {6 , ) , can be compared with observations. 



4-2 Determination of observed projected properties exhibiting correlations 

Determination of observed parameters which exhibit non-overlapping and well- 
separated correlations for different choices of intrinsic axial ratios is desirable 
problem since these parameters can be useful to get an estimate of intrinsic 
shapes of mass models, independent of viewing angles (SF94; TC01). It is also 
important to note here that correlation patterns between observed projected 
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properties evaluated at asymptotic radii, if exist, can save computing time due 
to availability of analytical expressions there. But, it would be meaningless to 
look for the correlations at asymptotic radii, since the observed properties at 
asymptotic radii are unobservable. With this point in mind, we computed ob- 
served projected properties at 0.25 and 4 effective radii, R e , in order to look 
for the above mentioned correlations at observable limits for the five differ- 
ent sets of (di, a,2, 03, 04) which include one as in Fig. 1 and remaining four 
as in Fig. 2. To compute R e , we fitted the azimuthally average surface den- 
sity of our models to i? 1//4 -profile ( de Vaucouleursl . 19481 ) between the interval 
0.25i? e < R < 4:R e and evaluated R e /b Q that depends on the intrinsic axial 
ratios and viewing angles for each of five adopted sets of constant parame- 
ters (01,02,03,04). In the spherical limit of our model, R e /b Q is found to be 
~ 4.4. However, it is also found that R e /b Q of our traixial models do not vary 
much from its value calculated in the spherical limit. Even then, we com- 
puted R e /b D for entire permitted range of intrinsic axial ratios and viewing 
angles for each set of (01,02,03,04) so that its exact value could be incorpo- 
rated in calculations of correlations. Furthermore, we defined new parameter 
(b/a)diff = (b/a).25R e — (0/0)4^, which exhibits ellipticity variation, and then 
found that the parameters {b/a)^R e and {b/a)diff are the better choice to rep- 
resent the desired correlations. 



In Fig. 4, we have shown the correlations between (6/a) 4 ^ e and (b/ajdifj for 
different choice of intrinsic axial ratios p = p Q = and q = q = Qoo roughly 
covering the parameter space of 1 > p > q > 0.5 when our model with con- 
stant parameters (ai, 02, 03, a 4 ) as in Fig. 1 is projected with varying viewing 
angles. We found that our model exhibits strong correlations which occupy 
well-separated regions in the parameters space of (6/a) 4j R e and (b/a)diff for 
different choices of (p, q) . 

It is interesting to note here that our model with remaining four values of 
(01,02,03,04) also show the qualitatively similar correlations as those found 
in Fig. 4. However, instead of considering all choices of (p, q) as in Fig. 4, we 
have presented them for only two choices of (p, q) = (0.65, 0.60) and (p, q) = 
(0.95, 0.60) in Figs. 5 and 6, respectively. The former choice of (p, q) represents 
highly prolate triaxial model, whereas the later choice of (p, q) corresponds to 
highly oblate triaxial one. Thus, these two choices of (p, q) are used to present 
non-overlapping correlations between a highly prolate and a highly oblate 
triaxial models that can be seen by comparing each plot of Fig. 5 with that 
of Fig. 6. 

Besides above found correlations where Ap = \p Q — p^l = and Aq = \q Q — 
goo I = are taken, it has also been seen that the correlation patterns still 
maintain as long as Ap and Ag are considered as 0.05. 

Finally, we have also examined the correlation plots by taking position angle 
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difference pa^f = (0* o — @*oo) as one of the observable parameters. Correlation 
plots between position angle difference paaif and (b/a)dif for a highly prolate 
triaxial and a highly oblate triaxial are presented in Fig. 7 where constant 
parameters (01,02,03,04) are same as in Fig. 1. Although some differences in 
the patterns between a highly prolate and a highly oblate triaxial are exhibited, 
they occupy overlapping regions. This suggests that pa di j is not a suitable 
observed parameter to develop well-separated correlations for two different 
choices of (p, q) required for the shape determination. 

But the non-overlapping correlations between the observed parameters (0/0)4^ 
and (b/a)diff are appeared to be useful in the prediction of the intrinsic shapes 
of the mass models by mean of the graphical method of shape estimate of SF94. 

Despite the radial profiles of projected properties of our models are much 
different from those of CTOO, its correlations between (6/0)4^ and (b/a)diff 
for different choices of (p, q) are found to be qualitatively similar to those 
shown in Fig. 1 of TC01 for CTOO model. This gives clear indication that 
our models can be good candidate for generating the ensembles of models, by 
considering different choices of Ap and Aq for each set of constant parameters 
(01,02,03,04), which could be useful in the technique of TC01 in order to 
get the model independent shape estimate of a model galaxy. To get more 
insight over this issue, an experiment to recover the intrinsic shape of a model 
galaxy is presented in the next section where we have used Bayesian statistics 
and exploited the correlations by employing the above mentioned observed 
parameters and ensembles of models. 



4-3 Recovery of the intrinsic shape of a model galaxy 



In or der to get an es timate of the intrinsic shape of a model galaxy, we fol- 
lowed Statle r (1994a||3, hereafter, ST94a,b) who used Bayesian statistics and 



considered photometric and kinematical data for estimating intrinsic shape of 
an individual galaxy. Moreover, he considered models of the type p = p(m 2 ), 
which give rise to concentric and coaligned ellipsoids, and used a constant 
value of ellipticity for a galaxy under consideration. 

In our present study of the shape estimate of a model galaxy, we do not use 
any kinematical data. Rather, we have incorporated the variations in ellipticity 
by considering (b/a)^ and (b/a)diff as the observed parameters which are 
calculated from our present model by using following intrinsic parameters: p^, 
goo, Ap, Aq, (9 ,<p) and (ai, a 2 , o 3 , a 4 ). Here and g^ are regarded as the 
shape parameters, while except of (6 ,<fi), the remaining intrinsic parameters 
are useful to generate ensembles of models. 

We launched a test case and considered an oblate traixial galaxy from our mod- 
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els with Pqo = 0.95, = 0.60, Ap = 0, Ag = 0, (#', 0') = (35°, 10°) and con- 
stant parameters (ai, a 2 , a 3 , a 4 ) = (0.2,0.9,0.2,0.9). The observed parameters 
of this selected model galaxy are found to be (o/a) /tp„ = 0.9398 and (b/a)dif f = 
—0.0597. We allowed a typical error of 0.01 (e.g.. lde Carvalho et al.l . ll99ll ) in 
ellipticities evaluated at 0.25i? e and at 4R e which in turns set errors of 0.01 
and 0.014 for (b/a)^ and (b/a)diff, respectively. 

In order to estimate the shape of this galaxy, we considered our present models 
with five different sets of constant parameters (ax, a>2, a^, a 4 ) as adopted in the 
preceding sections. For each choice of (01,02,03,04), we choose Ap > and 
Ag > in such a way that for each (p^, g^), (p Q , q Q ) are taken as (apoo, Pqoo) 
where choices of (a, 13) are (1,1), (1.05,1), (0.95,1), (1,1.05), (1,0.95). This 
gives the ensembles of models for a chosen value of {poo, goo)- For each model 
with a particular value of (poo, g^) and above given observed parameters and 
errors, we estimated the likelihood L of obtaining the observed data from the 
model and the posterior density P, which is the product of likelihood L and 
the prior density F, by the formulas given in ST94a,b. We assumed maximum 
ignorance of F = (sin 8 /4tt)F s by considering F s =constant (cf. ST94b). We 
then integrated the posterior density P over the viewing angles (6 , <ft ) and 
added it for different choices of (a, (3) and constant parameters (ai, a 2 , 03, a 4 ). 
We plotted the resultant posterior density as a function of (poo,goo) which is 
shown by dark-grey shade in Fig. 8. Here the darker shade corresponds to the 
higher posterior density. The inner and outer constant P contours, enclosing 
68% and 95% of the resultant posterior density, show lcr and 2cr error-bar 
regions, respectively. It is clear that the true shape of a model galaxy, in- 
dicated by a white cross, lies within the la region. We also found that the 
resultant posterior density is sharply peaked function of the shape parameters 
(Poo, <?oo) since the la region is sufficiently narrow. This satisfies the necessary 
condition of Bayesian statistics that the resultant posterior density to be rela- 
tively insensitive to the prior density in order to get the 'likelihood dominated' 
posterior density. 



5 Summary and discussion 

We have constructed a family of triaxial models by inclusions of additional 
radial functions, each multiplied by a low-order spherical harmonic, in the 
models of CT00. This procedure of building triaxial models by addition of 
the additional radial functions was initially suggested by TC01 to enlarge the 
ensembles of models in their technique of intrinsic shapes determination of 
model galaxies. Projected surface density £ of our present triaxial models can 
be calculated analytically which allow us to derive the analytical expressions of 
b/a and 0* at asymptotic radii. However, the approximate elliptical isodensity 
contours of S show the variations in b/a and G* as functions of R. This 
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suggests that our triaxial models can be considered as simple analytical models 
exhibiting ellipticity variations and isophotal twists. We have also studied 
the effects of the additional radial functions on the profiles of b/a and 0* 
by considering different choices of constant parameters (01,02,03,04). In this 
regard, different trends of the profiles of b/a and 0* are found for different 
contributions of additional radial functions. Furthermore, the profiles of b/a 
and 0* show small scale variations at inner R which are also seen in many real 
elliptical galaxies. Thus, our models can be useful to compare with observed 
photometric properties of real elliptical galaxies. 

TC01 pointed out that more ensembles of models, showing correlations be- 
tween the observed parameters for different choice of intrinsic parameters, are 
necessary to be explored to make their method of shape estimation more re- 
liable. In this regard, our models seem to be good candidate for providing 
another set of ensembles of models in the method of TC01, since result of 
shape estimation of a model galaxy using ensembles of our models is found 
to be satisfactory by utilizing the strong correlations between (b/a)±R e and 
(b/ a )diff- But it would be wise to look for larger sets of ensembles of models 
representing correlations between observed properties that could be useful in 
satisfactory determination of the shape of model galaxies before going ahead 
to apply this method of shape estimation to real ellipticals. 
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12 3 4 

R 



Fig. 1. Profiles of b/a, 0* and E as functions of R for a model with the in- 
trinsic parameters (p,q,9 ,<f>) = (0.65, 0.60, 60°, 10°) and the constant parameters 
(a\, a,2, 03, 04) = (0.2, 0.9, 0.2, 0.9). Here 0* is in degree, whereas R and E are in the 
units b D and M/Airb 2 , respectively. Like R, the constant parameters (01,02,03,04) 
are also given in the units of b Q . 
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R 

Fig. 2. Profiles of b/a as functions of R for four different choices of the constant 
parameters (01,02,03,04). The adopted values of (01,02,03,04) are given in top 
right corner of each panel, whereas the remaining intrinsic parameters of model, 
i.e., (p,q,9 ,</>), are same as those in Fig 1. Here R and constant parameters 
(ai, 02, 03, 04) are given in the units of b a . 
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Fig. 3. Profile of position angle 0* as a function of R. The intrinsic and constant 
parameters of model are same as those adopted in the bottom right panel of Fig 2. 
Units of R, 0* and constant parameters (01,02,03,04) are same as those given in 
Fig. 1. 
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Fig. 4. Each panel shows correlation between (6/0)4^ and (b/a)diff, when a model 
with a particular value of (p, q) is projected at various viewing angles. Considered 
values of (p, q) are indicated in top left corner of each panel, whereas constant 
parameters (01,02,03,04) are same as those adopted in Fig. 1. 
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Fig. 5. Correlation between (b/a)^ and (&/ a )di// for the fixed values 
(PiQ) = (0.65,0.60) but for different choices of constant parameters (a\, 02,03, 
as indicated in top left corner of each panel. 
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Fig. 6. Correlation plots as in Fig. 5 but for (p, q) = (0.95,0.60). 
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Fig. 7. Correlation plots between {b/a)dif and pa^if- The left panel represents cor- 
relation for model with (p, q) = (0.65,0.60), whereas the right panel shows it for 
(p,l) = (0.95,0.60). Constant parameters (01,02,03,04) are same as those adopted 
in Fig. 1. 
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Fig. 8. Shape estimate of a model galaxy using ensembles of models. A white cross 
represents the true shape of a model galaxy. Inner and outer contours represent la 
and 2a error-bar regions, respectively. 
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